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Abstract 

\ The inverse amplitude method is used to unitarize the two loop irir scattering 

' amplitudes of SU(2) Chiral Perturbation Theory in the I = 0, J = 0, I = 1, J = 1 

Q H > and I = 2, J = channels. An error analysis in terms of the low energy one-loop 

parameters ^1,2,3,4, and existing experimental data is undertaken. A comparison to 
standard resonance saturation values for the two loop coefficients 61,2,3,4,5,6 is also 
carried out. Crossing violations are quantified and the convergence of the expansion 
is discussed. 
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1 Introduction 

a: 

The 7T7T reaction is one of the theoretically cleanest processes in hadronic physics. This 
is due to the fact that crossing, unitarity and analyticity impose severe constraints to 
the scattering amplitude Q. Thus, a lot of attention has been paid to the study of this 
process and the determination of partial wave phase shifts P|-||10||. The current theoretical 
setup for such an approach is Chiral Perturbation Theory (ChPT) which is an Effective 
Field Theory embodying all these constraints, and leads to a perturbative expansion of 
the scattering partial wave amplitude in the (I, J) isospin-spin channel^ 

tiAs)=tf}(s)+tfj(s)+tfj(s) + --- (1) 

Here the expansion parameter turns out to be A 2 = m 2 r /(47r/ 7r ) 2 ~ 0.01 with m, = 
139.57MeV the physical pion mass and f w = 92.3MeV the weak pion decay constant. 
Hence tj) turns out to be proportional to A 2n . In the nir scattering case, the chiral 
expansion generates a hierarchy of corrections which depend on an increasing number of 
dimensionless and renormalization scale independent parameters. To lowest non-trivial 
order (LO,tree level) current algebra unique predictions for scattering amplitudes, tf)(s) 



in terms of f w and m r are generated [PL]. At next-to- leading order (NLO,one loop) four 



^-mail: jmnieves@ugr.es 
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3 We use the normalization such that the partial wave cross section is crij = (2J + l)(4ir/s)\tu(s)\ 2 
See also Eq. fl§) below. 
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k 


r(k,k) 


Set Ic (K14 decays) 


0.3 ± 1.2 


4.77 ±0.45 


2.9 ± 2.4 


4.4 ±0.3 


-0.69 


Set II (D— waves) 


-0.8 ±4.8 


4.5 ± 1.1 


2.9 ±2.4 


4.4 ±0.3 


-0.75 


Set III (Roy sum rules) 


-0.9 ± 1.2 


4.34 ±0.25 


2.9 ±2.4 


4.4 ±0.3 


-0.22 



Table 1: One loop h,2,s,i low energy parameters in ChPT for the parameter sets Ic, II and III of 
Ref. p|| as well as the correlation coefficient between T\ and fa. In the present paper we use only Set Ic 
and Set III. 



low energy parameters /1, 2,3,4 determine the amplitude t/j(s) [0, [13] • At next-to-next-to- 



leading order (NNLO, two loops) the amplitude tfj(s) can be expressed in terms of six 



parameters &i, 2,3,4,5,6 0> El- Unlike the one loop parameters /1, 2,3,4 which can be fixed 
from ChPT calculations confronted to experimental data from several sources, the two 
loop coefficients 61,2,3,4,5,6 are a bit more difficult to fix directly from experiment since 
the amount of data close enough to threshold is scarce. Because of this problem and 
motivated by the success of the resonance saturation hypothesis at the one-loop level and 
at a renormalization scale \x = 750 ± 250MeV |16|| , values for the two- loop 6's have been 



suggested on the basis of this hypothesis at that scale [jnj. This makes the, in principle, 
renormalization scheme independent fe's to have a spurious scale dependence. Since it 
is not really known what uncertainty should be assigned to this hypothesis, it has been 
suggested to ascribe a 100% error on the contributions to the low energy parameters 
determined by resonances 



m 



The numerical consideration of errors from ChPT requires taking into account consis- 
tent sets of low energy parameters, both ^,2,3,4 and 61,2,3,4,5,6, deduced from several sources 
18|, O, EH, |2l| . Moreover, there seems to be strong ant i- correlations between l\ and I2 on 



the light of two loop K u analysis [|T^, |2"0fl . This point has been thoroughly discussed in a 
previous work by two of us J^J and refer to it for further details. There, error propagation 
has been undertaken a la Monte Carlo, instead of using parametric statistics, by gener- 
ating a synthetic set of primary data. The basic assumption is that primary quantities, 
i.e., those obtained either directly from experiment or from an acceptable \ 2 fit, i.e., with 
a xVd.o.f. ~ 1, are Gauss distributed although perhaps with correlations. In practice, 
the distributions are represented by a population of finite but sufficiently large number 
(N = 10 4 ) of samples. The discussion in that work amounts to have three compatible 
sets of one loop I and two loop b low energy parameter distributions deduced from sev- 
eral sources. They are sumarized in Tables [I] and |2|. For consistency, we keep the same 
notation of our previous work. Set Ic in Ref. [^3|] corresponds to decays, following 
|T9"| and some statistical modeling designed to reproduce the fragmentary 

Set II corresponds to use D— waves as 



the fits of Ref. 
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information given in the analysis of Ref. 
proposed in the two loop 7T7t scattering calculation in standard ChPT |15l . Finally, Set 
III denotes the values of the low energy parameters obtained through Roy sum rules fol- 



lowing the lines of Ref. JT7 



231 that the anti-correlations between 



We have found in Ref. 

li and I2 persist in Set III, although they are not that strong. In the present work we 
discard Sets la and lb for being somewhat unrealistic. We also disregard Set II because 
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it produces too large errors as compared to Sets Ic and III. 



Despite of its great success, ChPT does not incorporate exact unitarity to a given 
order of the expansion and hence cannot account for resonances, in particular for the p 
and a resonances which appear in the I = 1, J = 1 and I = 0, J — channels respectively. 
To deal with this problem, several unitarization methods have been devised in the past 
in the context of tttt scattering ||24||-||36|| (for a recent and short review see e.g. Ref. P7j). 



In those methods, unitarity is restored in the different partial wave amplitudes, while 
crossing is violated [[ft], 3IJ. In the complex energy plane this corresponds to exactly 
take into account the unitarity right hand elastic cut ( s > 4m^) but to approximate the 
left hand cut ( s < ). Detailed quantitative studies reveal that the approximation used 
to take into account the left hand cut does not become critical to describe phase shifts in 
the scattering region, but it may significantly influence the violation of crossing and the 
values of the low energy parameters. Thus, there is some confidence that unitarization 
methods can indeed be used to enlarge the domain of applicability of ChPT to the study 
of intermediate energy hadronic reactions. Among these unitarization approaches the 
Inverse Amplitude Method (IAM) has successfully been applied to the description of 
meson-meson scattering incorporating up to one loop perturbative constraints. Original 
applications of IAM involved dispersion relation arguments [p7[ |29fl , which became rather 
cumbersome when incorporating coupled channels such as KK in tttt scattering. An 
algebraic derivation was soon found |32|] to provide an almost trivial generalization to the 



coupled channel case and a complete one loop analysis for all meson-meson channels has 
been carried out very recently in Ref. |33[. In addition to its very simple implementation 
from the standard chiral expansion, what makes this method particularly attractive is 
the fact that no new constants arise besides those already present in standard ChPT. 
Furthermore, the IAM method offers the possibility of systematic improvement according 
to the chiral expansion. Given the great success of this unitarization method in several 
meson-meson reactions including up to one loop corrections, there seems almost obvious 
to extend the calculation to the, in principle, more accurate description up to two loops. 
As we will show below, such an extension is not as trivial as one might think. Actually, 



there have been a previous calculation |3T| where an analysis of the IAM on the light 



of ChPT to two loops has been undertaken. The conclusion was that the IAM is a well 
converging scheme. It is fair to say that no effort was done to assign uncertainties in the 
low energy parameters, making somewhat hard to decide not only on the convergence 
itself but also on the compatibility with standard ChPT. 

Recently, theoretical restrictions for the s— wave scattering lengths have been obtained 
from an analysis of Roy equations E1J]. Unprecedented accuracy is obtained if in addition 
to the relativistic, crossing and unitarity demands from local quantum field theory, chi- 
ral symmetry constraints and the corresponding chiral expansion are implemented. The 
recent work on matching the Roy equation analysis |41| to the two loop ChPT ex- 



pansion fllj, [15[ has produced, using parametric statistics, the smallest error estimates 
for the low energy parameters, so far. Roy equations provide an extremely elegant frame- 
work to incorporate both crossing, analyticity and unitarity constraints on tttt scattering 
amplitudes. The set of non-linear inhomogeneous integral equations are not autonomous 
but require some high energy tails obtained from experiment as input. In the low energy 
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regime these so called driving terms can be described as polynomials which coefficients 
can be mapped to the low energy parameters of ChPT in the common region of valid- 
ity of Roy equations and ChPT. As a theoretical tool, Roy equations cannot be bitten 
by unitarization methods since the latter violate crossing to some extent. On the other 
hand, Roy equations have not been generalized yet to other processes different from tttt 
scattering and require a knowledge of high energy data which may not always be avail- 
able or accurate enough 0. Taking this fact into account and the time elapsed since the 



original work || and the recent update ETJ, it would be desirable but it seems unlikely 



that Roy equation techniques will become overnight a daily tool in hadronic physics. In 
contrast, unitarization methods based on ChPT require in principle no more work than 
ChPT itself which works well in the threshold region, but are able to describe, in addition, 
resonance physics and have been successfully applied to a variety of problems. Because of 
this the unitarization of tttt scattering amplitudes a la IAM provides a model case where 
we can learn about the virtues and drawbacks of the method and also on its convergence 
properties. 

In the present work we study the IAM of unitarization of the two loop ChPT ampli- 
tudes including a detailed error analysis based on the presently available information on 
the low energy parameters obtained from ChPT. By pursuing such a calculation we want 
to answer the question of whether or not low energy information plus unitarization repro- 
duces the data beyond the domain of applicability of standard ChPT. In common with 
other unitarized calculations it is not clear how to avoid the unavoidable and prejudiced 
choice of a particular unitarization method. Moreover, given the unitarization method it 
is hard to estimate uncertainties due to higher orders in the expansion. Our only hint so 
far, is to compare successive orders in the scattering phase shifts with their corresponding 
error-bars and determine whether or not practical convergence requirements are met. We 
do this analysis using Set Ic and Set III of one- and two loop low energy parameters 
£12 34 and 61,2,3,4,5,6 respectively of our previous work p3|. 



The paper is structured as follows. In Sect. [] we provide some basic definitions in order 
to fix notation. We also estimate unitarity violations and the failure of a perturbative 
definition of phase-shifts to describe the data in the region above threshold. In Sect. |3| we 
analyze the IAM phase-shift predictions as well as the corresponding threshold parameters, 
together with estimates on the amount of crossing violations in terms of Roskies sum 



rules |38|, [39] inherent to any unitarization method. Motivated by previous experiences 
such crossing violations can be amended by a suitable generalization of the IAM. This 
point is analyzed in Sect. |3.5| . Although, there has been some work on one loop IAM of 
unitarization, we present in Sect. [| an updated analysis from the point of view of the 
convergence of the expansion for the unitarized phase shifts. Finally, we draw our main 
conclusions in Sect. [5[ In the Appendix |S| we provide some information not presented in 
our previous paper |23[ and also relevant for the present work such as correlation matrices 



of both low energy constants and threshold parameters. 



See however the recent work on irN scatteringBj 
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2 ChPT to two loops and unitarity violations 



2.1 Basic definitions 

Let tjj(s) be the partial wave scattering amplitude for the reaction tttt — > tttt at the 
Center of Mass (CM) energy y/s in the I J isospin-spin channel 

e 2iSu(s) _ j 



with a(s) = Jl — Ami/ s the CM momentum and d~u(s) the corresponding phase-shifts. 
Two particle unitarity corresponds to real 5u(s) and can be written as a non-linear 
relation on the amplitude 

Imtuis) = a(s)\t u (s)\ 2 (3) 

or, equivalently, as a linear relation on the inverse amplitude 

lmtj}(s) = -a(s) (4) 

For a tttt scattering amplitude calculated in the chiral expansion sketched in Eq. ([!]), the 
lowest order amplitude tf)(s) is a real function for S— and P— waves, and vanishes for 
D— and higher ones. The NLO amplitude ijj(s) develops an imaginary part for S— and 
P— waves but becomes real for D— waves, and so on. The exact unitarity relation of 
Eq. requires, at a perturbative level, the set of relations 

Imtfj(s) = (5) 
Im4 4 j( S ) = a(s)\tf)(s)\ 2 (6) 
lmtf)(s) = 2a(s)tf)(s)Retf)(s) (7) 

Standard ChPT fulfills exact crossing symmetry at any order of the expansion, but violates 
unitarity. Two aspects are related to this violation. In the first place, a necessary condition 
for unitarity is the fulfillment of the inequality 

a(s)\t u (s)\ = | sin 5 u (s)\ < 1 (8) 

This unitarity limit yields values of s a bit too highP]. A better way to quantify the 
unitarity violation is by defining the quantity 

\Ju(s) = \l + 2ia(s)t IJ (s)\ (9) 

Below the two pion production threshold \fs = Am^ ~ 560MeV the elastic unitarity 
condition requires Uu(s) = \e 2l6lJ ^\ = 1. Strictly speaking, if unitarity is violated, 
Uu(s) 7^ 1, there is no way besides perturbation theory for a real phase shift to fulfill 

5 For instance, for the isoscalar S— wave one gets at LO the inequality fulfilled in the range s < 
^y/irfir ~ 670MeV. We will show below that unitarity violations take place at significantly lower energies. 
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simultaneously Eq. (P and Eq. (0). Expanding Eq. (|) according to Eq. ([[]) the standard 
ChPT phase shift may be computed yielding 



a(2), 



+ (T{8) 



tfks) 



2ia(s)tf)(s)tf)(s) 



+ ... 



(10) 



The elastic unitarity condition corresponds to Sfj PT (s) being real, which is automatically 
fulfilled if the perturbative unitarity relations, Eq. (|7|), are used, and one effectively gets 



^ hPT (*) 



°( s )t?j{s) + a{s)Retfj{s) + a(s 



( 4 ), 



Retf)(s) + la(s)H?)(sr 



+ ...(11) 



Close to threshold, the scattering amplitude can be written in terms of the threshold 
parameters, scattering lengths, a/j, and slopes, bu, defined by 



tu(s) = 2m n (s/4 - ml) J [a u + b u (s/A - m*) + ■ ■ 



(12) 



The scattering amplitudes tjj(s) present kinematical zeros of order J at s = 4m^.. Chiral 
symmetry implies the existence of dynamical zeros for the S— waves, named chiral or 
Adler zeros. In ChPT, Adler zeros may be determined perturbatively, i.e. 



tlj(s A ) = tf)(8 A ) + tT>(s A ) + tfj(s A ) + 



f(6)/ 



.(2) 



,(4) 



,(6) 



Expanding the solution s A = s A + s A ' + s A ' + 

-- 



we getf] 



s (4) 



■■ -4 4 i(^ 2) )/4 2 i(4 2) )' 

,(4)/ (2)y,(4)/ (2)s 
(6) _ ,(6)/ (2)w (2)/ (2)y tlj {S A ) tjj {S A ) 
A — L IJ\ b A )l L IJ\ b A ) i" r (9) , r21x.no 



[4 2 j( S ; 2) vi 2 



At lowest order, non-kinematical zeros are located at 



(13) 



(14) 
(15) 

(16) 



J2) 
S A 

J2) 
S A 



-ml 1 = J = 
2ml 1 = 2 J = 



(17) 
(18) 



From the formulas above, direct application of ChPT requires a scrupulous separation 
between different chiral orders. As we have already mentioned above, the NLO amplitudes 
tjj(s) depend linearly on four dimensionless parameters /i, 2,3,4- These parameters are 
supposed to be independent of /„- and m n and therefore they are zeroth order in the chiral 

For a general function there would be a term involving the second derivative of t\]. This term 

(2) 

disappears from the formulas since t\j is a linear function of s. 
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Figure 1: Unitarity condition for standard NNLO ChPT amplitudes in tttt scattering for 
S— and P— waves defined by Ujj(s) = \l+2ia(s)t IJ (s)\ . Upper panel: Set Ic of Ref. [fffi]. 
Lower panel: Set III of Ref. [^] . In the calculation the parameters 6° defined in Eq. (|19|) 
and given in Table § have been used. The unitarity condition requires Uu(s) = 1. 

expansionQ Finally, at NNLO the amplitude can be expressed in terms of six independent 
parameters &i,2,3,4,5,6- As it has been shown in the original two loop calculation works 



14], |j| these b— coefficients contain a zeroth order piece and a second order piece 

k = 6° + Ah t (19) 

This makes the separation a bit subtle because, when plugged into Eq. ([!]), a unwanted 
eighth order correction is induced. From the point of view of ChPT this contribution has 
to be dropped out as far as the complete eighth order calculation is not available. On the 
other hand, these higher order corrections are numerically small as can be deduced from 
Table 0. 



7 Practical calculations require, however, a truncation of the chiral expansion and confrontation to ex- 
perimental data, and hence some higher order systematic uncertainties remain in the one loop parameters 
besides the experimental uncertainties 
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s" 2 (MeV) 



s" 2 (MeV) 



Figure 2: Standard NNLO ChPT phase shifts (in degrees) for tttt scattering for S— and 
P— waves after Eq. ([II]). Upper panel: Set Ic of Ref. Lower panel: Set III of 

Ref. [p3 |. In the calculation the parameters b® defined in Eq. (pT9|) and given in Table |2] 
have been used. Combined data from Refs. [S]-[10|. 



2.2 Numerical results 



The CM energy dependent figures with error-bars presented in this work are generated 
as follows. If we have an energy dependent function, F, in terms of a set of random 
parameters, (ai . . . a n ), distributed according to some statistical law, a random variable 
for any fixed value of s, F(s] ai, . . . , a n ) is generated. Obviously, for a non-linear parameter 
dependent function the mean value of the curve is not equal to the function of the mean 
values, (F(s )) 7^ F(s; (ai), . . . , (a n )). There is nothing wrong with this and one 

could simply bin the distributions for any fixed s value. Nevertheless, to make the results a 
bit more portable we wish to quote such a function of the mean values, F(s; (ai), . . . , (a n )) 
as our central curves. To assign an upper and lower error bar (the distribution may in 
general be asymmetric) relative to this central value, we bin the distribution and firstly 
exclude the 16% top values and the 16% bottom values of the distribution. The remaining 
bins comprise the 68% of the distribution values, the distance from the upper and lower 
values to our central value provide the upper and lower error-bars respectively. Evidently, 
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Figure 3: Unitarity condition for standard NNLO ChPT amplitudes in tttt scattering for 
S— and P— waves defined by Ujj(s) = \l + 2ia(s)tjj(s)\. Upper panel: Set Ic of Ref. 



Lower panel: Set III of Ref. |23[ . In the calculation the parameters 6j defined in Eq. (|IJ 
and given in Table § have been used. The unitarity condition requires Uu(s) = 1. 



our bands correspond to a 68% confidence level. This procedure of assigning errors fails 
for extremely asymmetric distributions, such that the central value turns out to be within 
the discarded upper or lower 16% intervals. Although we find that this situation seldom 
takes place in our calculation, in such a case we proceed in a different way. We first 
discard the 16% upper and lower intervals and then compute the arithmetic mean, which 
we assign as the central value. To control on the quality of this second definition of central 
value we always compute, whenever possible, both definitions and find that the differences 
are numerically less significant than the error-bars. 

The results for the unitarity condition, Uu(s) as defined in Eq. (||), can be seen in 
Figs. [I] and |3] in terms of the 6° and the hi = b® + Abi coefficients respectively and for 
the parameter Sets Ic and III. As we see, standard ChPT theory violates unitarity in a 
systematic manner well below the resonance region including uncertainties for Set Ic. For 
Set III only the P— wave exhibits this behaviour, due to large uncertainties in the S— wave 
unitarity violation. The perturbatively defined phase-shifts, Sfj PT (s) as given by Eq. ([IT]) 
can be seen in Figs. 0and|4] in terms of the 6° and the 6j = b^ + Abi coefficients respectively 



9 




s" 2 (MeV) 



s" 2 (MeV) 



Figure 4: Standard NNLO ChPT phase shifts (in degrees) for tttt scattering for S— and 
P— waves after Eq. flIH|). Upper panel: Set Ic of Ref. P3fl. Lower panel: Set III of 



Ref. [p3 |. In the calculation the parameters hi defined in Eq. ( |I~9D and given in Table 
have been used. Combined data from Refs. [S]-[10|. 



and for the parameter Sets Ic and III. As we see from the figures, the perturbatively 
defined phase shifts seem compatible with experimental data whenever the corresponding 
unitarity condition is compatible within uncertainties with elastic unitarity. Let us remind 
that threshold parameters for Set III have uncertainties similar or a bit smaller than Set 
Ic [p3] . This also holds for higher energies but there appear some systematic discrepancies 



with the data, slightly favouring Set III. 

Threshold parameters for Set Ic and Set III at the two loop level have been computed 
in our previous work J23[. There, a separation of the tree- level, one loop and two loop 



contributions with their corresponding error estimates have been presented. 

The partial S— and P— wave amplitudes in the unphisical region below threshold 
and above the left cut, < < 4m^, are depicted in Fig. |!| In this region partial wave 
amplitudes are real and present real single zeros. The single zero at 4m^ in the P— wave 
is of kinematical origin as can be seen from Eq. (p!2|). However, zeros in the S— waves 
are dynamical consequences of chiral symmetry. As we see, the agreement between the 
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b° 


10 

°5 
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Set Ic 


-9.1 ±2 


8.2 ±1.7 


0.3 ±0.3 


0.66 ±0.07 


c 7+3.2 
°- ' -3.9 
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Set III 


-10.71IJ 


9.8 ±1.7 


-0.16 ±0.40 


0.58 ±0.04 
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2 +0 - 8 
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Ab 4 


Ab 5 


Ah 


Set Ic 


-2.4+gj 


3.0 ±0.3 


-o.5i8:? 


0.19 ±0.04 








Set III 


-2-4185 


o «+0.4 
z -°-0.3 


-0.4 ±0.2 


1 ^+ - 15 

u - io -0.04 









Table 2: Two loop 61,2,3,4,5,6 low energy parameters in ChPT for the parameter Sets Ic and III 
of Ref. |23| ] used in the present paper. The relation to the 61,2,3,4,5,6 parameters used in the work is 
bi = l&Trbi. We also show explicitly the decomposition bi — 6° + Abi referred to in Eq. (|TjJ). 



SA/ml 


(I, J) = (0,0) 


(J, J) = (2,0) 


ChPT-Ic Fig 


•1 


0.38(6) 


2.03(5) 


ChPT-III Fig 


I 


0.43(6) 


2.00(5) 


ChPT'-Ic Fig 


•1 


0.38(6) 


2.03(5) 


ChPT'-III Fig 


•I 


0.43(6) 


2.00(5) 



Table 3: Non-kinematical Adler zeros for 5— wave I — and 1 = 2 amplitudes for ChPT and the 
parameter sets Ic and III of Ref. |23) . We also indicate the phase-shift figures which correspond to these 
zeros. Errors are given in brackets. 

parameter Sets Ic and III is very good within uncertainties. The two-loop location of 
Adler zeros with error estimates is given in Table |3|. The additive structure of the ChPT 
amplitude makes a numerically small distinction between making the separation bi = 
b^ ± Abi. As we see, the isotensor S— wave chiral zero does not move within uncertainties 
from its tree level value of Eq. ( |T8"D both for Set Ic and Set III. For the parameter Set III, 
the two-loop shift of the S— wave isoscalar Adler zero this is almost the case. Nevertheless, 
for parameter Set Ic there is a systematic shift of about 20%. 

3 IAM to two loops and perturbative matching 
3.1 Algebraic derivation 

The idea of the method is quite simple and we review it here for the sake of completeness. 
If instead of expanding the amplitude, one considers the inverse amplitude and expands 
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(l,J)=(0,0) 








, 1 























s/m, 




Figure 5: S 1 — and P— partial wave amplitudes fu(s) = \fstu(s) (in fm) for nir scattering 
in standard ChPT to two loops in the unphysical region, < s < 4m^. Upper panel: Set 
Ic of Ref. |p3[j . Lower panel: Set III of Ref. ||23|| . Normalization is such that at s = 4m^ 
one has for the S— waves the scattering length. In the scattering region ( s > 4m^. ) this 
figure corresponds to Fig. 0. 



according to the chiral expansion (assuming tjj(s) ^ ), one gets 



tu(s) 



a(s) cot 5 j j 



ia(s) 



4%) 



+ 



4 6 j( s ) 



+ ... (20) 



One may check that the unitarity relation, Eq. (|j), is exactly preserved due to the pertur- 
bative relations of Eq. (|7|). Note that a direct application of the IAM including up to two 
loops can only unitarize S— and P— waves. To unitarize D— waves a three loops calcula- 
tion would be needed. Since such a calculation has not yet been done, we will restrict to 
S— and P— waves in the present work. The structure of Eq. fl2"Dp makes possible to have 
poles in the second Riemann sheet, i.e. zeros of t _1 (s), but this is done at the expense 
of some fine-tuning between several orders. Actually, the IAM method assumes that the 
inverse amplitude, l/tu(s), is small, which is particularly true in the neighbourhood of a 
resonance. 
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3.2 IAM Phase-shifts 



The best way to quantify the goodness of a unitarization scheme such as the IAM is 
to check whether or not the information contained in the low energy parameters, in 
conjunction with the unitarized amplitude given by Eq. (p0|) , predicts within acceptable 
errors the phase shifts in the region above the threshold. To proceed further we have to 
fix in some way our sets of parameters. An alternative, and actually complementary point 
of view is to make a direct fit to the data. Unfortunately, this involves a 10-parameter 
fit, and moreover there are some parameters, like for instance 1$ and I4, for which tttt 



scattering is not very sensitive. Actually, as it has been recognized in Ref. there are 
two kinds of low energy parameters according to the properties of the corresponding terms 
in the partial wave amplitudes, 

• Class A Terms that survive in the chiral limit comprising l±, I2, r§ and r§ 

• Class B Symmetry breaking terms corresponding to the remaining low energy pa- 
rameters I3, Z4, ri, r 2 , r 3 and r 4 . 

where the parameters determine the pure two loop contribution to the amplitude. On 
the basis that chiral symmetry breaking is a small effect, we expect a higher sensitivity 
of the scattering data to variations of the class A parameters. 



3.2.1 Naive scheme 

The simplest and most direct way to look at the quantitative predictions of the IAM is 
to take the unitarized amplitude Eq. fl20| ) for all partial waves and propagate the errors 
in the one- and two-loop parameters £1,2,3,4 and 61,2,3,4,5,6 respectively. As we have already 
pointed out before, this may be a dangerous procedure since the two loop parameters 
contain a higher order piece, but one might argue that since the numerical effect on the 
6's should be small (see Table ^), one might expect an overall small effect anyhow. Along 
these lines, we present in Fig. |6| the results obtained by using the parameter Sets Ic and 
III of Ref. |[23|| . As can be seen from the figures, the errors are huge and there is even a 
trend to discrepancy in the p channel for Set III. 



3.2.2 Monte-Carlo scheme 



Having one realized of the dangers of making a naive Monte Carlo error propagation from 
the analysis of Ref. ||23|| , we proceed now in a different manner. We consider Sets Ic 
and III of Ref. |23| for both the one loop /^-parameters and the zeroth order two loop 



parameters 6° as explicitly given in the Appendix of Ref. ]T5|]. The results are shown in 
Fig. [?]. Although in the S— waves there are no big differences as compared to Fig. || 
in the p channel, the effect at intermediate energies makes the predicted phase-shifts not 
only compatible with data but also the error-bars get substantially reduced. 
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Figure 6: IAM Unitarized phase shifts (in degrees) for nir scattering for S— and P— 
waves. Naive scheme (see main text). Upper panel: Set Ic of Ref. Lower panel: Set 
III of Ref. pi. Combined data from Refs. 0-pL 



3.2.3 Partial Fit scheme 



The results of the two previous schemes suggest that some fine tuning mechanism is 
needed, since the large difference between them is due to either considering the b param- 
eters as a whole or to explicitly split them as different orders, according to bi = 6° + Afej. 
This fact seem to indicate that their numerical values may be constrained to a large extent 
by performing a x 2— fit. As we have mentioned above this involves 10-parameters. An- 
other possibility can be displayed by making selective fits in some subsets of parameters 
propagating the errors in the remaining ones. After trying out several combinations, we 
have found that, indeed, the low energy parameters of class A, i.e., those not vanishing 
in the chiral limit, are enough for a satisfactory fit to the data. In practice, the proce- 
dure is as follows. For either Set Ic and Set III, we generate a sufficiently large sample 
(N = 10 4 proves large enough), of class B parameters, i.e., those corresponding to chiral 
symmetry breaking terms in the scattering amplitudes. For any member of the class 
B parameter population a x 2 ~ fit °f class A parameters is performed. This procedure 
yields distributions for li, I2, and parameters whence straightforward error analysis 
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Figure 7: IAM Unitarized phase shifts (in degrees) for tttt scattering for S— and P— 
waves. Monte-Carlo scheme (see main text). Upper panel: Set Ic of Ref. [2~3|. Lower 



panel: Set III of Ref. [23. Combined data from Refs. M-M 



may be undertaken. The choice of these parameters rather than the 6's has the advantage 
that the separation b,i = 5° + A6j may be explicitly done within the fitting procedure, 
and hence a correct bookkeeping of chiral orders in the inverse amplitude is implemented. 
The result of the fit isfl 

h = -0.141°$ f 2 = 44±0.1 10 4 r 5 = 1.07l[J;i 10 4 r 6 = -0.35l°J| (21) 

which produces x 2 /d.o.f. = 69.9/ (67 — 4) = 1.11, a rather satisfactory value as can also be 
clearly seen in Fig. |[ Besides, we obtain the correlation coefficient r(l%, h) = 0.22, There 
are no Set Ic and Set III labels because l\ and 1%, which provide this label, are determined 
from the fit. Both class A and class B parameters contribute to the total errors. The 
errors corresponding to fitted (class A) low energy parameters can be determined by 
employing the standard procedure of changing the \ 2 from its minimal value by one unit. 
We find that they are rather small, so the quoted errors in Eq. (|21|) are dominated by 



8 We apply the following energy cuts and scattering data. For the isoscalar S— wave we cut at y/s — 
610MeV the data of Refs. [§, |, |, |, @, §. For the isotensor S- wave we cut at 970MeV the data of 
Refs. H, ||. For the isovector P— wave we cut at MeV. 
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Figure 8: IAM Unitarized phase shifts (in degrees) for tttt scattering for S— and P— 
waves. Partial fit scheme (see main text). There are no Set Ic and Set III labels because 
li and l 2 , which provide this label, are determined from the fit. Combined data from 
Refs. Uncertainties in the curves stem from those of class B parameters and the 

scale for which resonance saturation is assumed, /i = 750 ± 250MeV. 



the uncertainties in the class B low energy parameters and the scale at which resonance 
saturation is assumed, \i = 750 ± 250MeV. 

As we see, the out-coming values of the fitted parameters, particularly l\ and T 2 , 
are very much in agreement with the standard ChPT estimates of Ref. [fffi| . It is also 



interesting to note that the values of r§ and are consistent within error-bars with those 
assumed from resonance saturation provided with a 100% uncertainty as suggested in 



Ref. []l7fl • Transforming these values into b parameters we get 



6i = -12.1±i;l b 2 = ll.Btli h = -0.2918:1? 

6 4 = 0.75 ±0.02 b 5 = 3.14l°;g 6 6 = 0.55±g;|| (22) 

These numbers have been constructed by adding the out-coming 6° to the second order 
Abi contribution, although both contributions enter the fit in a non-additive way. As we 
see, the resulting values agree with those expected from standard ChPT analyses within 
uncertainties, with the sole exception of b$ which turns out to be inconsistent. The reason 
is that the corresponding value for r 6 comes out exactly opposite sign as that expected in 
resonance saturation. 

Finally, we present in Fig. [| the IAM unitarized partial S— and P— wave amplitudes 
in the unphysical region < s < 4m^. using the Monte Carlo scheme for both Sets Ic 
and III. Obviously, the kinematical zero of the P— wave at s = 4m^. remains fixed. On 
the other hand, the non-kinematical Adler zeros of the S— wave amplitudes do not move 
from their lowest order locations given by Eq. flT8] ) but become higher order zeros, as 
can be observed from the figures in the small bumps around the zeros and analytically 
in Eq. (|20"D. Although this effect is undesirable, we see that from a direct comparison 
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Figure 9: S— and P— partial wave amplitudes fu(s) = y/stu(s) (in fm ) for mr scattering 
in the IAM to two loops in the unphysical region, < s < Am\. Monte-Carlo scheme. 
Upper panel: Set Ic of Ref . [^3| . Lower panel: Set III of Ref . [23] . Normalization is such 
that at s = Am\ one has for the S— waves the scattering length. In the scattering region 
( s > Am\ ) this figure corresponds to Fig. |7|. 



of the standard ChPT amplitudes of Fig. [5] with the IAM unitarized ones of Fig. |9] one 
may conclude that the violation of the order of the non-kinematical zero does not have 
quantitative dramatic consequences. 



3.3 Crossing violations 

As we have said, the IAM restores unitarity but violates crossing symmetry. The in- 
teresting point is that although unitarization, which has to do with the right-hand cut, 
may provide a satisfactory energy dependence and, hence a description of the data in the 
scattering region, it can still do so with low energy constants (LEC's) which differ from 
those expected in standard ChPT. The precise numerical values of the LEC's depend on 
how the left-hand cut is handled and approximated before or after unitarization. On the 
other hand a proper left-hand cut is a direct consequence of crossing symmetry in the 
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s, t, u representation. It is therefore interesting to study these crossing violations on a 
quantitative basis. A particular way of doing this at the level of partial wave amplitudes 
is by using Roskies sum rules |38|, 39 1. Two ways at least have been introduced in the 
literature to characterize crossing violations in a quantitative way. Defining the quantities 
introduced in Ref. |30] 



c 2 
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These relations can be written in the general form 
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Crossing symmmetry implies C\ = 0. We also consider the definitions introduced in 
Ref. H 
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Table 4: Roskies sum rules violations in percentage as defined by Eq. ( p7| ) and introduced in Ref. pQ] 
for the IAM method and the parameter sets Ic and III of Ref. ^T|. We also indicate the phase-shift 
figures which generate these violations 
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Crossing symmetry implies in this case A4 — B\ = for i from 1 to 6. Formally, if the 
unitarized amplitudes embody ChPT to some order, these sum rules will be identically 
verified to the same order. In previous works, the numerical fulfillment of the sum rules 
has been tested to one and two loops, but no error estimates have been taken into account. 
Due to this, some ad hoc dimensionless crossing violations have been defined 0, [HJ. In 
Refs. [§t], 0] the ratio 



FL = 100 x 



Jo"" ds EijUu^tujs) 
So m * ds \J2u wjj,i(s)t/j(s)| 



(26) 



is introduced^ whereas in Ref. ||0| the violation 

A: - Bt 



V, = 100 x 



Ai + Bt 



(27) 



is defined. The advantage of providing error bars to the crossing sum rules C{ = or 
A4 — Bi = is obvious; the dimensionless quantity is naturally defined as the size of 
the uncertainty relative to the mean value. This allows to make a definite statement on 
crossing violations in terms of statistical uncertainties. Nevertheless, we quote in Tables f| 
and H the Vi and Ri values of Ref. [|(| and Refs. |3(J |34|] respectively since they give an 
idea of how large are these violations in percentage terms. 



As can be inferred from Table |], crossing violations as introduced in Ref. do not 
seem to be dramatically large, although this depends on their particular definition. The 
most serious violations appear in the V5 rule, which combines both isospin 5— wave chan- 
nels and the P— wave. The computed uncertainties provide a less pessimistic impression, 
since in some cases these violations are compatible with zero. Moreover, if the crossing 

9 We use specifically the formula proposed in Ref. since it is unambiguous. The formula of Ref. 
is misleading, probably due to a misprint. 
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Table 5: Sum rules violations in percentage as defined by Eq. ( p6| ) and introduced in Refs. |30|, Qfor 
the IAM method and the parameter sets Ic and III of Ref. |23| . We also indicate the phase-shift figures 
which generate these violations 





a 00 m v 


booml 


10 • anml 


10 • b n ml 


10 • a 2 0^7r 


10 • &20«^ 


ChPT-Ic Fig. | 
ChPT-III Fig. | 


0.214(5) 
0.208(6) 


0.27(1) 
0.25(1) 


0.37(1) 
0.374(8) 


0.06(1) 
0.053(7) 


-0.42(1) 
-0.44(1) 


-0.76(2) 
-0.80(2) 


ChPT'-Ic Fig. | 
ChPT-III Fig. | 


0.214(5) 
0.208(6) 


0.27(1) 
0.25(1) 


0.37(1) 
0.374(8) 


0.06(1) 
0.053(7) 


-0.42(1) 
-0.44(1) 


-0.76(2) 
-0.80(2) 


IAM-Ic Fig. | 
IAM-III Fig. | 


0.220(7) 
0.211(8) 


0.30(2) 
0.27(2) 


0.37(8) 
0.37(6) 


0.048(6) 
0.046(5) 


-0.42(1) 
-0.44(1) 


-0.76(2) 
-0.80(2) 


IAM'-Ic Fig. 
IAM'-III Fig. 


0.221(8) 
0.213(8) 


0.29(2) 
0.27(2) 


0.38(1) 
0.381(9) 


0.072(1) 
0.064(1) 


-0.42(1) 
-0.44(1) 


-0.76(2) 
-0.80(2) 


FIT-III Fig. | 


0.216(5) 


0.280(7) 


0.376(6) 


0.058(5) 


-0.43(1) 


-0.79(1) 



Table 6: Scattering lengths, 
parameter Sets Ic and III of Ref. 1 23 
threshold parameters. 



ai.j_an& slopes bu defined by Eq. (|T^) for the IAM method and the 
We also indicate the phase-shift figures which correspond to these 



violation definition of Refs. |3(| [34| is evaluated we see from Table |5] that these sum rules 
are better fulfilled. The effect of uncertainties on the violations has been overlooked in 
previous works |3(| 31]. Nevertheless, we point out that generally speaking there are 
systematic, though small, crossing violations. In the partial fit scheme, corresponding to 
Fig. |] the crossing violation defined in Refs. 
than 0.1%. 



30, p4| are compatible with being smaller 



3.4 Threshold parameters 

The chiral expansion is expected to work best at low energies. But even so, threshold 
parameters such as scattering lengths ajj and effective ranges bjj defined by Eq. (|T2|) turn 
out to get corrections at each order of the expansion. The IAM method is constructed to 
reproduce ChPT at all energies but in the lowest orders of the l// 2 expansion. Thus, if 
we go to the threshold region we do not exactly reproduce the standard ChPT behaviour. 
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Table 7: Roskies sum rules violations in percentage as defined by Eq. ( p7j ) and introduced in Ref. 
for the generalized IAM method of Ref. p3| and the parameter sets Ic and III of Ref. |p3| . 



Generalized IAM 


Rx 


R2 


R 3 


R4 


R 5 


Set Ic Naive 
Set III Naive 
Set Ic Non-Naive 
Set III Non-Naive 


0.0818:8! 
0.0818:81 
-o.n±8:8l 
-0.0818:81 


-0.06+8-1 

0.02+8! 
-o.3l8':l 
-o.3l8:l 


-o.o4l8:?8 
-o.o5±8:8! 

o ns +0 - 10 
n n4+ - 10 


-o.o6l8:8i 

-0.0318:81 
_n 1 ?+ - 03 

u - iz -0.04 

n 1 n +003 


-0.0518:8^ 
-o.oo4±°:8I 

n 9+0.05 

u - z -0.06 

-0.218:? 



Table 8: Sum rules violations in percentage as defined by Eq. ( p6|) and introduced in Refs. 
the generalized IAM method of Ref. Q and the parameter sets Ic and III of Ref. [E3I . 



for 



Nevertheless, as can be seen from the figures the difference of the standard ChPT ampli- 
tudes and the unitarized ones is actually very small. Our results for the two loop IAM 
threshold parameters are presented in Table || The fact that ChPT and ChPT' entries 
of the table are the same within errors is not accidental; it merely reflects the additive 
combination hi = 6° + Abi and the smallness of Ab{. A more detailed table containing 
the explicit separation in tree-level, one-loop and two loops contributions as well as 68 % 
ellipses of the S— wave scattering lengths can be found in Ref. 031 . In general we see that 



the IAM threshold parameters are compatible within errors with the ChPT ones. The 
only exception is the slope bu in the p channel for the two Monte-Carlo schemes. The 
partial fit scheme provides compatible an and bu values with slightly better accuracy 



3.5 Generalized IAM 



Roskies sum rules provide a set of necessary conditions for a crossing symmetric tttt 
scattering amplitude. It is has been noticed that the IAM method transforms the non- 
kinematical single zeros of the partial wave amplitudes into + 1-order zeros of the IAM 
unitarized amplitudes, being N the order of the chiral expansion (see the denominators 
[tfj{s)] N+1 of Eq. (|20|) ) . Since the integrals involve the interval < s < Arn^ between 
the right- and left-hand cuts, this higher order zeros clearly influence the fulfillment of 
the crossing sum rules. However, there is no unique way how to modify the chiral zeros 
behaviour in order to achieve a better fulfilment of crossing. To overcome this difficulty 
several interesting methods have been proposed effectively rectifying the amplitude be- 
haviour in the unphysical region, though one should say that none of them is entirely 
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satisfactory from the point of view of the mathematical properties that one wants to im- 
pose a priori on the amplitude. In Ref. |3(| it has been suggested ( scheme II of that 
work) to use a dispersion relation for the inverse amplitude, tu(s)~ l . In such a way, 
not only the unitarity cut but also the position of the single chiral zero, which becomes 
a single pole for the inverse amplitude may be enforced from the beginning. The left- 
hand cut is assumed to be that one of ChPT up to a certain negative s = —A 2 value ( 
A 2 = 0.5 — 0.6GeV 2 ) and a constant up to s = — oo. This procedure has the disadvantage 
of introducing a new variable (the cut-off A) into the problem not present in the ChPT 
amplitude. In addition, it does not take the shift in the non-kinematical chiral zeros into 



account, and imposes the tree level ones. More recently, in Ref. [[34[] a generalized I AM 



(2) 

to two loops has been proposed. If so is the chiral Adler zero to lowest order t\j [sq j 
then the following expression for the inverse amplitude is suggested pl | 
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This expression violates exact unitarity since 
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and, in addition, has a single zero at the lowest order approximation of the chiral zero. 
The slope coincides with the one obtained in ChPT as can be seen from the formula 



tu(s) 
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(30) 



In the limit 4j ( s o)+4j ( s o) — ¥ in Eq. (pS|) the generalized IAM of Ref. |34| reduces to the 



standard IAM method of Eq. (EO) and also unitarity is exactly fulfilled. In practice, both 
unitarity violations and the absence of a shift for non-kinematical zeros are numerically 
small. It has been shown in Ref. |34[ that the generalized IAM improves the fulfillment 
of the Roskies sum rules, but no uncertainties estimates were considered. Using the two 
definitions of crossing violations given by Eqs. (|26|) and (j27|) suggested in Refs. |3(| |34| 
and [5C] we show in Tables |7] and ^ respectively that this is indeed the case, provided 



uncertainties in the parameters are taken into account. 



4 On the convergence of the IAM method 

The IAM method can be systematically implemented to any order in the chiral expansion 
with no additional LEC's than those required by standard ChPT. There arises the question 
to what extent is this method convergent. To answer this question in practice we can 
only compare one-loop and two loop predictions for the unitarized phase shifts. Such 
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Figure 10: Unitarity condition for standard NLO-ChPT amplitudes in irir scattering for 
S— and P— waves defined by Ujj(s) = \l+2ia(s)tu(s)\ . Upper panel: Set Ic of Ref. [f2^ 
Lower panel: Set III of Ref. p3) The unitarity condition requires Ujj(s) = 1. 



a comparison makes sense only if errors in the LECS are also transported, as we have 
repeatedly done along this work. Actually, in Ref. [|29| the one loop error analysis of the 
IAM phase shifts was estimated by varying the low energy constants. In this section we 
reanalyze this question by using the updated values of the one loop coefficients given by Set 
Ic and Set III taking into account by means of a Monte-Carlo simulation the important 
anticorrelations between l\ and I2 determined in Ref. |23|| . Following the same systematics 
of the two loop calculation, we show in Fig. [K] the unitarity condition of Eq. (Q). As one 
would expect, unitarity violations of one- loop ChPT occurr at lower energies. The NLO 
ChPT phase shifts, defined through Eq. ( |10|) are depicted in Fig. |ll|- The general trend 
follows a similar pattern to the two loop calculation, although some important differences 
emerge. Firstly, the uncertainties in the phase shifts are smaller at NNLO than at NLO 
in the threshold region, as one would expect from the fact that threshold parameters are 
more accurately determined at NNLO than at NLO in the region above threshold 

10 This circumstance is not trivial and it only happens for Setslc and III. The Set II of Ref. |^3| shows 
cases where predictive in the threshold parameters is lost, and the NNLO result is no more accurate than 
the NLO. 
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Figure 11: Standard NLO-ChPT phase shifts (in degrees) for tctc scattering for S— and 
P— waves after Eq. flIH|). Upper panel: Set Ic of Ref. P3| . Lower panel: Set III of 
Ref. |3[. Combined data from Refs. @- 13 . 



the situation is exactly the opposite, the two loop calculation produces larger uncertainties 
than the one loop one. In addition, by comparison of Fig. [11] and Figs. |2| and (|4] in the p 
channel) one sees that the discrepances in the region above threshold are larger than the 
estimated uncertainties, with an overall trend to improvement in the two loop calculation. 
A similar trend, although in a less clear manner, is observed in the two S— waves. The 
unitarity condition, Eq. (g), gives us a good idea on the applicability of standard ChPT to 
one and two loop approximations. Nevertheless, the agreement of the perturbative phase 
shifts, Eq. (|TUD, with experiment seems to extend up to a region where the unitarity 
violation may be as large as 10 — 20%. 

The one loop IAM phase shifts are depicted in Fig. [12|. Apparently, the general picture 
provided by NLO-ChPT looks better than what one obtains by comparing any of the two 
Monte-Carlo schemes studied in Sect. |3] based in NNLO-ChPT. By looking at any of 
these two-loop schemes, Fig. ^| and Fig. [7], we realize that there is a clear loss of predictive 
power; the errors in the two loop phase shifts are larger than the discrepancy between 
their mean value and the one loop mean value. Finally, in Fig. [IB] a partial one-loop fit 
procedure in l\ and I2 parameters to the data is presented, where variations on the [3 and 
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Figure 12: NLO-IAM Unitarized phase shifts (in degrees) for tttt scattering for S— and 
P— waves. Monte Carlo scheme (see main text). Upper panel: Set Ic of Ref. [ESI. Lower 
panel: Set III of Ref. J23|. Combined data from Refs. @-|g. 



U are taken into account. The result of the fit is 



-0.44 ± 0.02 



U = 5.51 ±0.04 



-0.81 



(31) 



where the errors reflect the uncertainites in Z3 and Z4. Here x 2 /d.o.f. = 191/(67 — 2) = 
2.94 almost three times larger than in the two loops case (1.11), and too large to be 
considered a satisfactory drescription of the scattering data. Such a large x 2 value makes 
the determination of the uncertainties of l\ and I2 due to the error bars in the fitted data 
meaningless. As we see, the obtained values for the fitted parameters are compatible 
with the corresponding two loop partial fit procedure, Eq. ( pi]) , although the errors in 
the one-loop case due to uncertainties in Z3 and Z4 are much smaller than in the two-loop 
case. This may be an indirect consequence of the large \ 2 value. 
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Figure 13: NLO-IAM Unitarized phase shifts (in degrees) for tttt scattering for S— and 
P— waves. Partial fit scheme (see main text). The errors in the curves are due to 



uncertainties in Z 3 and Z 4 . Combined data from Refs. p|j-fT0|l 

5 Conclusions 



In the present work we have presented a thorough study of the Inverse Amplitude Method 
to unitarize the NLO (one-loop) and NNLO (two-loop) ChPT nn scattering amplitudes 
below the KK threshold. To this end, we have considered several one-loop /i, 2,3,4 and 
two loop &i, 2,3,4,5,6 parameter sets along the lines discussed in our previous work . 
Particularly interesting in this work is the role played by the uncertainties in these pa- 
rameters. To complement the analysis and provide some quantitative motivation we have 
determined unitarity violations within standard ChPT, with error estimates. They take 
place at much lower energies than the unitarity limit suggests. Moreover, we have also 
shown the systematic discrepancy with the data in the region above threshold if phase 
shifts are defined perturbatively. The discussion is complicated by the fact that the two 
loop parameters 61,2,3,4,5,6 ma Y by splitted in a zeroth order contribution 6? 2 3 4 5 6 an d a 
higher order contribution A6 12 , 3,4,5,6 which slightly spoils the chiral counting. The dis- 
tinction caused in the phase shifts by including or not the higher contribution is small 
within ChPT. Motivated by this we have unitarized the two loop amplitude, and devised 
several schemes to predict the phase-shifts from threshold up to the resonance region. 
The effect of consistently treating or not A6 12 ,3,4,5,6 as higher order is much stronger for 
the IAM unitarized phase shifts. Typically, a factor of two difference or larger in the 
uncertainties is encountered. In any case, they are rather large, although seem consistent 
with the scattering data. This indicates a kind of fine tuning going on and suggests a 
fit to the data to determine the low energy parameters which remain in the chiral limit, 
keeping the remaining low energy parameters within their error bars. The result of the 
fit is satisfactory, although there appears a discrepancy in the 66 coefficient. Neverthe- 
less, the predicted partially fitted phase shifts vary within very small uncertainties, not 



far from the recent ChPT analysis of Roy equations carried out in Ref. El. Despite of 



these features, the IAM produces crossing violations which have been quantified in terms 
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of Roskies sume rules. Generally speaking, they are not very large in percentage terms, 
although in some cases the uncertainties are so large that no conclusion can be drawn. 
We have also studied some proposals to generalize the IAM in order to achieve a better 
fulfillment of crossing properties. Finally, we have adressed the convergence of an expan- 
sion based on the IAM and increasing order of approximation in ChPT. By comparing 
NLO (one-loop) and NNLO (two-loop) IAM predicted phase shifts we see at the present 
stage a lack of predictive power; the errors in the two loop phase-shift are larger than the 
difference between the central one-loop and two-loop values. This is a direct consequence 
of the low accuracy in the two loop parameters. 



A Correlations among NNLO low energy constants 
and threshold parameters in standard ChPT 

The correlation matrix, defined as usual, 

r ij = ( X i X j) 

Cj - (Cj) 

1 N 

(f(Ci,...,C n )) = —Y.f( C ^,---,C n ,a) (32) 

iV (2=1 

being Cj any of the low energy constants or threshold parameters is provided below as 
obtained from our N = 10 4 finite size samples both for Set Ic and Set III. Taking into 
account the central values and their errors given in Table |2| and Table |^ and ignoring 
the slight error asymmetries the parameter Sets are fully portable by going through di- 
agonalization to principal axis in parameter space and making a Monte-Carlo Gaussian 
simulation in each principal direction. 

Set Ic 



r(bi, bj 



/+1.00 










\ 


-0.74 


+1.00 










+0.53 


-0.77 


+1.00 








-0.49 


+0.64 


-0.49 


+ 1.00 






+0.09 


+0.04 


-0.41 


+0.10 


+1.00 




V-0.08 


+0.23 


-0.54 


+0.23 


+0.57 


+ 1.00/ 


&00 


a u 


a 20 


boo 


bn 


b 2 o 
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r Thrcshold 

ij 



/+1.00 










\ 


+0.13 


+1.00 










+0.58 


-0.21 


+1.00 








+0.75 


-0.13 


+0.37 


+1.00 






+0.02 


+0.89 


-0.07 


-0.10 


+1.00 




V+0.20 


+0.06 


+0.57 


+0.38 


+0.11 


+ 1.00/ 



(34) 



Set HI 



r{b h bj) 



^Threshold 

ij 



/+1.00 










\ 


-0.73 


+1.00 










+0.53 


-0.76 


+ 1.00 








-0.12 


+0.20 


-0.14 


+ 1.00 






+0.22 


-0.16 


-0.19 


+0.21 


+1.00 




V+0.06 


+0.02 


-0.29 


+0.19 


+0.47 


+ 1.00/ 




a u 




»oo 


hi 




/+1.00 










\ 


-0.08 


+1.00 










+0.67 


-0.46 


+ 1.00 








+0.81 


-0.32 


+0.55 


+ 1.00 






-0.11 


+0.82 


-0.23 


-0.18 


+ 1.00 




V+0.37 


-0.45 


+0.71 


+0.59 


-0.21 


+ 1.00/ 



(35) 



(36) 
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